High order harmonic generation in noble gases using plasmonic field enhancement 
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We present theoretical studies of liigh-order harmonic generation (HHG) in rare gases driven 
by plasmonic field enhancement. This kind of fields appears when plasmonic nanostructures are 
illuminated by an intense few-cycle laser and have a particular spatial dependency, depending on 
the geometrical shape of the nanostructure. We demonstrate that the strong nonhomogeneous 
character of the laser enhanced field plays an important role in the HHG process and significantly 
extends the harmonic cutoff. Our models are based on numerical solution of the time dependent 
Schrodinger equation (TDSE) and supported by classical and semiclassical calculations. 
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I. INTRODUCTION 

One of the most prominent demonstrations of the non- 
hnear character of the interaction of a strong laser with 
matter is high-order harmonic generation (HHG) [l|, Q . 
In the past two decades, this phenomenon has been an 
important subject in both theoretical and experimental 
physics, as a reliable tool for generating coherent light 
sources in the ultraviolet (UV) to extreme ultraviolet 
(XUV) spectral range. HHG based coherent light sources 
open the possibility for extracting temporal and spatial 
information with attosecond and sub- Angstrom resolu- 
tion i. As a results, HHG has found numerous appli- 
cations in basic research, biology, material science, and 
even lithography [1]. 

The physical mechanism behind HHG can be under- 
stood by a semi-classical three-step model [ll, |B| • Accord- 
ing to this model, in the first step, an electron reaches 
the continuum with no kinetic energy through tunneling 
ionization. In the second step, the released electron sub- 
sequently accelerated in an oscillating electric field until 
it is driven back by the laser field toward its parent ion 
or molecule. In the third step, upon its return, it recom- 
bincs with the core and emits a high-energy photon, with 
energy equal to the sum of the electron kinetic energy and 
the ionization potential. In this model, the propagation 
is treated classically while the ionization and recombina- 
tion processes are described quantum mechanically. The 
maximum photon energy (cutoff) producible with high 
harmonic generation is given by Ip + 3.17C/p, where Ip is 
the atomic ionization potential, while Up ~ I/iuj'^ is the 
ponderomotive energy in the laser field of intensity / and 
frequency ui. 

For high-order harmonic generation to occur, one needs 
a laser field with intensity greater than 10^'^ W/crn^), two 
orders of magnitude larger than the output of the current 
modest femtosecond oscillator. It means an additional 
process like chirped-pulse amplification is needed to reach 
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the required intensity for generation of high harmonics 
using noble gases. Even so, the XUV based on HHG has 
low duty cycle and efficiency 

The recent demonstration based on surface plasmon 
resonances as light amplifiers could overcome these dif- 
ficulties 0, Hi- By maneuvering surface plasmon reso- 
nances, the laser electric fields can locally be enhanced 
by more than 20 dB g ^ without the need of extra 
cavities or laser pumping process. Indeed, in such a sys- 
tem the intensity of the locally enhanced field become 
large enough to exceed the threshold intensity for har- 
monic generation in noble gases. In particular, in the 
case of gold bow-tie-shaped nanostructures, the initial 
laser field (800 nm laser with intensity 10^^ W/cm^) en- 
hanced sufficiently to generate the 7th (114 nm) to the 
17st (47 nm) harmonics. Moreover, each nanostructure 
could potentially act as a point-like source to generate 
high harmonic radiation. This would allow coUimation of 
these coherent radiations by means of constructive inter- 
ference. This would open a wide range of possibilities to 
spatially rearrange nanostructures to enhance and even 
shape the HHG spectra 0| • Additionally, one of the most 
unwanted effects of the harmonic propagation, namely 
the phase matching, appears to be not important at all 
in the HHG using plasmonic fields Q. 

The physical mechanism of HHG based on plasmon- 
ics can be explained as follows (for full explanation see 
01 )■ A femtosecond low intensity laser pulse couples to 
the plasmon mode and initiates a collective oscillation 
among the free charges within the metal. This causes a 
large resonant enhancement of the local field inside and 
at the nanostructure vicinity. This enhancement is well 
above the threshold intensity for generating high har- 
monics. Consequently, by injecting rare gases into the 
site of the enhanced field HHG can be produced. In 
here, the enhanced field is spatially inhomogeneous in 
the region where electron dynamics take place. In addi- 
tion, the movement of the electron in the enhanced field 
is restricted in space. These two features imply strong 
modifications in the harmonic spectra, as was shown by 
several authors Furthermore, the influence of 

the nonhomogeneous character of the laser electric field 
in the photoelectron spectra was recently assessed [l^ . 
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The outcome of Kim ct al. Q experiment, however, has 
been recently under intense scrutiny |16l4l8| . Between 
the points to elucidate are the real intensity enhance- 
ment of the input laser field and the damage threshold of 
the gold bow-tie nanostructures. According to the finite- 
element simulations presented in Ref. 0, the local field 
intensity could be enhanced by 4 orders of magnitude (40 
dB). It means for input laser intensities of the order of 
10^^ W cm~^ the enhanced field intensity could reach to 
the order of 10^^ W cm~^ in the vicinity of the bow-tie 
nanostructure. On the other hand, from the measured 
high-order harmonic spectra it appears that these num- 
bers are unrealistic. The highest measured harmonic was 
n « 17, a limit corresponds to intensities of the order of 
5 X 10^'^ W cm~^ (we have estimated this value using the 
well established three step or simple man's model 0). 
It is clear that the measured value for enhanced field 
intensity is different from the ones obtained in the finite- 
element simulations of Ref. 0] ■ In our theoretical models 
we have considered this point using a reduction factor 
in the field enhancement obtained from the finite ele- 
ment simulations. Moreover, for this work we use a laser 
pulse with a longer wavelength, in order to have larger 
electron excursion, to explore the effect of the field non- 
homogeneity more clearly. 

To our best knowledge, apart from Ref. Q there is only 
one more experiment in which HHG was produced from 
noble gas using metallic nanostructures Q. Recently, 
however, different kinds of nanostructures, e.g. nanopar- 
ticles or nanotips, have been used to generate high energy 
electrons directly from the nanostructure in the absence 
of the gas and to study the distinct and new character- 
istics of these nanosytems (see e.g. [T9l - [23j ) 

The paper is organized as follows. In the next section 
we present a numerical calculation of HHG by incorpo- 
rating the actual functional form of the electric field in 
the vicinity of a metal bow-tic shaped nanostructure in 
the time dependent Schrodinger Equation. A classical 
and semiclassical analysis of an approximated model for 
nonhomogencous fields is presented in Section 3. We fin- 
ish with our conclusions and a brief outlook. 



II. SINGLE ATOM RESPONSE OF PLASMONIC 
ENHANCED FIELDS 

Most of the numerical and semiclassical approaches 
to study laser-matter processes in atoms and molecules, 
in particular high-order harmonic generation (HHG), 
are mainly based on the dipolc approximation in which 
the laser electric field (E(r, i)) and its vector poten- 
tial associated (A(r, t)) are spatially homogeneous in 
the region where the electron dynamics takes place, i.e. 
E(r,t) = E(t) and A(r,<) = A{t) On the other 

hand, the fields generated using plasmonics nanosystems, 
such as metallic nanoparticlcs, nanotips and gold bow- 
tie shape nanostructures, present a strong spatial de- 
pendence, which cannot be ignored. From a theoretical 



viewpoint, the HHG process using homogeneous fields 
can be tackled using different approaches and approxi- 
mations (for a summary see e.g. H^, and references 
therein). In this section we employ the Time Dependent 
Schrodinger Equation (TDSE) in order to study the har- 
monic radiation generated by a model atom when it is 
illuminated by a spatially inhomogcncous electric field. 
This in an actual field, which is obtained directly from 
3D finite-element (FE) simulations where the real pa- 
rameters of the bow-tie shaped nanostructures, as those 
used in 0], are considered. For linearly polarized fields, 
which is the case of our study, the dynamics of an atomic 
electron is mainly along the direction of the laser electric 
field and as a result it is a good approximation to cm- 
ploy the Schrodinger equation in one spatial dimension 
(ID-TDSE) [l| which can be written as follows: 



^ dt 



1 



(1) 



VL{x,t) 



where Va{x) is the atomic potential and VLix,t) repre- 
sents the potential due to the laser electric field. In here, 
we use for Vatom the "quasi- Coulomb" or "soft-core" po- 
tential 



Va{x) 



^x^ + e 



(2) 



which first was introduced in [26| and since then has been 
widely used in the reduced dimensions studies of laser- 
matter processes in atoms. The required ionization po- 
tential can be defined by varying the parameter ^ in Eq. 
([2]). The laser potential VL{x,t) of the laser electric field 
E{x,t) is given by 



VL{x,t) = E{x,t)x. 



(3) 



In Eq. ([3]), the spatial dependency of E{x,t) can be 
defined in terms of a perturbation to the dipole approxi- 
mation and it reads 



E{x,t) = Ea fit) (1 + h{x)) sinwt. 



(4) 



which is linearly polarized along the x-axis. In Eq. Q, 
i?o, w and f{t) are the peak amplitude, the frequency 
of the coherent electromagnetic radiation and the pulse 
envelope, respectively. In addition, h{x) represents the 
functional form of the nonhomogencous electric field and 
it can be written as a power series of the form h{x) = 
Tlii^i biX^- The coefficients bi are obtained by fitting the 
actual electric field that results from a FE simulation 
considering the real geometry of different nanostructures 
(See panel (b) of Fig. 1). In this work we use for the 
short laser pulse a trapezoidal envelope given by 



fit) = 



(t3-t2) 



for 0<t <ti 
for ti <t<t2 
for t2 <t <t3 
elsewhere 



(5) 



where ti = 27rno„/a;, ^2 = ^1+ 27rnp/a;, and ^3 = ^2 + 
2TTnoff/uj- rion, n-p and rioff are the number of cycles of 
turn on, plateau and turn off, respectively. 

We employ ^ = 1.18 in Eq. ^ such that the binding- 
energy of the ground state of the ID Hamiltonian coin- 
cides with the (negative) ionization potential of Ar, i.e. 
£gs = —15.7596 eV (—0.58 a.u.). Furthermore wc as- 
sume that the noble gas atom is in its initial state (ground 
state (GS)) before we turning the laser {t — —00) on. 
Equation ([T]) is solved numerically by using the Crank- 
Nicolson scheme [l|. In addition, to avoid spurious re- 
flections from the spatial boundaries, at each time step, 
the electron wave function is multiplied by a mask func- 
tion [23|. 

The harmonic yield of an atom or molecule is propor- 
tional, in the single active electron approximation, to the 
Fourier transform of the acceleration a{t) of the active 
electron H^. That is. 
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dte~'"*a(t) 
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with a{t) is obtained by using the following commutator 
relation 
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In here, 'H(t) and 'ii{x,t) arc the Hamiltonian and the 
electron wave function defined in Eq. ([T]), respectively. 
The function D{uj) is called the dipole spectrum, which 
gives the spectral profile measured in HHG experiments. 
For solving Eq. ([1]), the gap size g of the gold bow-tie 
nanostructure is taken into account restricting the spatial 
grid size (see Figure 1). 

The electric field intensity distribution inside the gap 
of gold bow-tic shaped nanoantcnnas was computed nu- 
merically bySD Finite Element Method (COMSOL Mul- 
tiphysics) I29l . using the gold optical properties taken 
from Ref. [30| . The nanoantenna is formed by two iden- 
tical (isosceles) triangular gold pads (longest altitude of 
600 nm and the smallest acute angle of 30°) separated 
by an air gap g. The apices at corners were rounded (10 
nm radius of curvature) to account for limitation of cur- 
rent fabrication techniques and avoid nonphysical fields 
enhancement due to tip-effect. The out of plane thick- 
ness is set to 25 nm. These parameters yield to a dipo- 
lar bonding resonance centered at around A = 1800 nm 
when considering gaps ranging between 12 nm and 15 
nm. This particular value of A was chosen according to 
the availability of laser sources [3l|. On the other hand, 
the selected laser wavelength allows the electron to have 
excursions of the order of the gap g and consequently to 
confine its movement. Classically the electron excursion 
in an oscillating electric field is given by the so-called 
quiver radius uq, which is (x ^fl\^ where / is the laser 
intensity. For instance, in here, for intensities / of ~ 10^^ 
W cm~^, ofQ can have a value about ±80 a.u. (±4.5 nm). 

The panel (a) of Fig. 1 displays the calculated electric 
field enhancement in the gap of the bow-tic structures 



FIG. 1: (Color online) Panel (a) electric field enhancement 
in the gap region of the gold bow-tie nanostructure extracted 
from the FE simulations; panel (b) functional form of the elec- 
tric field E{x, t) where the solid lines is the raw data obtained 
from the FE simulations and the red dotted line is a nonlin- 
ear fitting; panel (c) high-order harmonic generation (HHG) 
spectra for Ar with ionization potential Sgs ~ —0.58 a.u., 
laser wavelength A = 1800 nm and intensity I = 1.25 x 10^* 
W-cm~^ at the center of the gap 2; = 0. We use a trape- 
zoidal shaped pulse, Eq. ([S]), with Hon = 3, Uoff = 3 and 
Tip = 4 (about 60 fs). The gold bow-tie nanostructure has a 
gap g = 15 nm (283 a.u.). Black line indicates the homoge- 
neous case while magenta line indicates the nonhomogeneous 
case. The arrow indicates the cutoff predicted by the semi- 
classical model for the homogeneous case 



when illuminated by a linearly polarized (x-axis) plane 
wave at 1800 nm. The field-enhancement profile is ex- 
tracted for the bow-tie long axis through the middle of 
the gap, so the successive problem is reduced to ID as is 
shown in panel (b) of Fig. 1. Additionally we normalize 
the electric field by setting E{Q,t) = 1. We observe a 
typical amplification between 30-40 dB, i. e. 3-4 orders 
of magnitude between the input intensity and the inten- 
sity at center of the gap. In a real experiment, however, 
the enhancement can be smaller than our calculations; 
nevertheless, for a similar system, it was shown that one 
can obtain values of more than 20 dB 0, Q . 

Panel (c) of Figure 1 depicts the harmonic spectra for 
a bow-tie shaped nanostructure with gap g = 15 nm and 
for a laser wavelength of A = 1800 nm considering an 
homogeneous electric field, i.e E{x,t) = E{t) and a non- 
homogeneous electric field using Eq. (j?]). The enhanced 
laser intensity is / = 1.25 x 10^^ W cm~^ at the center of 
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the spot {x = 0). In order to reach this intensity value, we 
consider plasmonic field enhancements between 25 and 
35 dB with input intensities in the range of 4 x 10^^- 
4 X 10^° W cm~^. These intensities would be well below 
the damage threshold of the nanostructure employed (see 
e.g. 0)- In both cases, we use a trapezoidal shaped pulse 
with three optical cycles turn on (no„ = 3) and turn off 
(rioff = 3) and a plateau with 4 optical cycles {up ~ 4), 
i.e. 10 optical cycles in total which is about 60 fs. 

For the homogeneous case we have an harmonic cutoff 
at around 204ti; as shown by the arrow in panel (c) of 
Fig. 1. In fact, our calculation are in excellent agree- 
ment with the semiclassical model Q . For the nonhomo- 
geneous case, however, we observe a substantial increase 
in the harmonic cutoff, which is about 50 % higher than 
the cutoff generated by a homogeneous electric field. This 
new feature emerges due to the combination of the non- 
homogeneous character of the electric field and the con- 
finement of the electron motion [l^. In fact, the laser 
ionized electron could be absorbed by the metallic surface 
(for details see [ll|) and in this way only the short elec- 
tron trajectories would be allowed. This circumstance 
is comparable to the restriction in the electron motion 
incorporated in our theoretical model (32j . 



III. CLASSICAL AND SEMICLASSICAL 
APPROACHES 

In this section we present classical and semiclassical 
studies of a particular case of nonhomogeneous field. In 
particular, we only consider the linear term of the series 
for h{x) (see previous Section). This field represents, for 
instance, a first order approximation of the near-field in 
the surrounding of a metal nanoparticle [s^ [s^. The 
main advantages of this approximation is the possibility 
to perform a quasi analytical approach using the Strong 
Field Approximation (SFA) [61 . In this way, the main 
features behind the quantum trajectories are elucidated 
and characterized. 



as 

Eeff{t,x) E{t)il + ax), (9) 

where a <C 1 is a parameter that characterize the 
strength of the inhomogeneity. 

By substituting ^ into ([5]), we have 

x{t) = -E{t){l + 2ax{t)) = -E{t, x). (10) 

Classically, the electron trajectory can be found by 
solving Eq ([TU)) . In here, we solve this equation by ap- 
plying the Picard iteration [ssf method and restrict our- 
selves to the first order (for more details see [l^l)- In the 
SFA, the action is defined in terms of the vector poten- 
tial field A{t), the counterpart of the velocity x{t) in the 
classical equation, which corresponds to the integration 
of x{t) with respect to t. By assuming that the electron 
starts its movement from the origin with zero velocity, 
i.e. a;(0) = and i(0) = , we will have 

x{t) - A(t) - A(to) - A(to)(t - to). (11) 

with X{t) = J^^dt'A{t'). Furthermore, if we assume the 
potential field is zero at time to, then 

x{t) = J dt'A{t'). (12) 

By using Eqs. and ([TU)) . the effective vector potential 
along the electron trajectory Aeff(t,x) reads 

A,ff{t)^A{t)+2aAait), (13) 

where 

Aa(t)^J dt"A{t")-J dt"A^{t"). (14) 



A. Strong Field and the Newton equation 

We use Eq ([3]) in order to calculate the Newton equa- 
tion of motion of an electron in a nonhomogeneous elec- 
tric field. In particular, the Newton equation of motion 
for an electron in this field is given by 

i{t) = -\/,Vl = -x\/Aff{t,x) - E,ff{t,x), (8) 

where Eeff{t, x) gives the effective electric field along the 
electron trajectory. For the homogeneous case, in which 
the laser field does not have spatially dependency and it is 
just E{t), the Newton equation reads x{t) = —E{t). On 
the other hand, if the spatial dependence of the enhanced 
laser electric field is perturbativc and linear with respect 
to position, then the effective field can be approximated 



B. Transition Amplitude 

The strong field approximation (SFA) is based on the 
assumptions that the influence of the laser field is ne- 
glected when the electrons are bound to their target 
atoms and the binding ionic potential is neglected when 
the electrons are in the continuum. As a result, the free 
electrons in the continuum are described by field-dressed 
plane waves, which are known as Volkov states (36l . Wh . 

Based on the well established model of Lewenstein [y], 
the SFA transition amplitude for HHG reads (in atomic 
units) 

/OO f^t /■ 

dt dt' d^kd;,,{k{t))d,on{k{t')) 
-OO J — OO J 

g-iS(n,k,M') + c.c. (15) 
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with the action 



sin,k,t,t') = JJ^±Ail]Ldr + ip{t-t)-nt + 

2aJ^Aa{T)[k + AiT)]dT+ (16) 

2a^ J^ Al{T)dT 
and the prefactors 

d^on{Ht')) = {k{t')\H,nt{t')\<l>o) (17) 



(18) 



In here, Ip, fc, ft , Hint{t'), Odip and Cx denote the ioniza- 
tion potential of the field- free bound state |(/io), the drift 
momentum of the electron in the continuum, the har- 
monic frequency, the interaction of the system with the 
laser field, the dipole operator and the laser polarization 
vector, respectively. 

Eq. (|15p describes a physical process in which an elec- 
tron, initially in a bound state |0o) with energy /p, in- 
teracts with the laser field by Hint{t') at the time t' and 
ionized into a Volkov state |k(i)). Subsequently, it prop- 
agates in the continuum from time t' to t. At the time t 
it is driven back to the parents ion and recombines with 
the core emitting high-harmonic radiation of frequency 
51. The SFA is not gauge invariant, thus, Eqs (fTT]) and 
([T5|) will have different form in the length and velocity 
gauges. For the length gauge, the interaction Hamilto- 
nian and the Volkov wave function are given by -ff (i) = 
r.E{t',x) and k(t) = k + yle//(i), respectively, while for 
the velocity gauge they are Hl^^{t) = [k + Aeff{t')]/2 
and k(t) ~ k, respectively. In this work, we consider 
the length gauge and assume that the electric field is lin- 
early polarized along x-axis. Furthermore, we consider 
an hydrogenic Is state for the field-free bound state |^o)- 

As a result, Eqs. ^ and (0 yield 



(k(t')) « 



drec{^{t)) OC 



(19) 



(20) 



respectively. 



C. Saddle-point equations 



The transition amplitude (llSp was computed employ- 
ing the saddle-point method [Sg. For this approach 
one needs to obtain the saddle points for which the ac- 
tion is stationary, i.e. for which dtS{Q,^'k,t,t') = 



dfS{n,\i,t,t') = dkS{n,k,t,t') = 0. The solutions of 
the saddle point equations are directly related to the clas- 
sical trajectories, which allows to demonstrate the cutofi^ 
of the HHG. The stationary conditions upon t,t' and k 
lead to the saddle-point equations 



[k-f A(i')] -2e/3(t') = -2/p 



dT[k + A(r)] + 2ari{T) = 



and 



with 



and 



2al3{t) 



riir) 



Aa{T)dT. 



(21) 



(22) 



(23) 



(24) 



(25) 



Eq. ([21]) gives the conservation law of energy for the 
electron tunnel ionized at the time t' . Eq. ([^^ constrains 
the intermediate momentum of the electron as well as 
guarantees that the electron returns to its parent ion. 
Finally, (Eq. I23p expresses the energy conservation of 
the electron at the recombination time t and generation 
of a high frequency photon 57. 

The terms l3{t'), r]{t) and l3{t) in the Eqs. (EH, (EH), 
and (j23p . respectively, give the nonhomogeneous char- 
acter of the laser field and they are zero for the homo- 
geneous case (e = 0). For small e, the solutions of the 
saddle equations are generally complex since Eq. (ET|) has 
no real solutions, unless Ip — >■ 0. For nonhomogeneous 
field, 3.17[/p, where Up = E'^/(Auj'^) is the ponderomotive 
energy, does not directly corresponds to the maximum ki- 
netic energy that the electron gains in the continuum . 
Depends on the nonhomogeneous character of the field, 
for the positive and negative Aa this value is larger and 
smaller than 3.17C/p, respectively. 

We now examine the drift momentum k of the electron 
at the time of the tunneling . For that, we take the limit 
in which the electron reaches the continuum with zero 
kinematical momentum, /p — > 0. As a result, Eq. (PTj) 
reads 



-A{t') + eAa{t'){lTVi) 



(26) 



In here, k has two different solutions with one lowering 
and the other exceeding and the homogeneous drift mo- 
menta. This result is different from the homogeneous 
case, in which k = —Ait'). This deviation emerges from 
the strength of the inhomogcneity e and the shape of the 

Aa{t'). 
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D. Orbits 

111 this section, we investigate the contribution of indi- 
vidual trajectories to the HHG cutoff of the homogeneous 
and nonhomogeneous cases by examining the solutions 
of the saddle point equations. To get a better insight 
into the problem, we employ a monochromatic field with 
E{t) = E()Siii(ujt)ex, where is the polarization vector 
along the x axis. By using the conventional relationship 



A(t) = - f E(t')dt' 



(27) 



and applying some trigonometric identities, the electric 
field (Uni) and potential field ^ read 



m 



Eosm{ujt){l + 2afi'm{ojt)/uj'^), (28) 
AoCos{Lut) +2aAait), (29) 



respectively, where Aq ~ Eq/lu and 

Aa{t) ^ sin(wt)/4a; - J^^tj^. 



(30) 



In terms of the pondermotive energy of the homogeneous 
field, the drift momentum of Eq. ([^5]) yields 



k = -2y;7pCos(wt') + e — sin(wt') - 2?/pi' (1 T \/3) 

(31) 

Based on the above equations, we solve the saddle point 
equations defined in Eqs. ([?!]) - in terms of the ion- 
ization (i') and recombination (t) times. 
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FIG. 2: (Color online) Dependence of harmonic order on the 
release time i! and the recombination time t of the electron 
for the homogeneous field (e = 0). We consider an hydro- 
gen atom, for which the ground-state energy is /p = 0.5 
a.u., in a linearly polarized, monochromatic field of frequency 
u = 0.057 a.u. and intensity 7 = 3 x 10"W/cm^ Panels (a) 
and (b) give the harmonic order as a function of the ioniza- 
tion and recombination times of the SFA model, respectively, 
while panels (c) and (d) depict the the harmonic order in 
terms of ionization and recombination times of the classical 
calculations, respectively. The red dashed lines correspond to 
the harmonic cutoff. 



In this paper, we present these solutions briefly, for 
more detailed analysis we refer to our previous paper [s^ . 
In Figs. [2]and[3l we plot the harmonic order as function 
of the real parts of the ionization (t') and recombination 
(t) times for the case with e = and e = 0.003 (panels 
a and b for t' and t, respectively), respectively. In these 
figures, in comparison to the SFA model, we also present 
the classical solutions of t' and t (panels c and d, respec- 
tively) . From these figures it is clear that the SFA resem- 
bles the classical calculations. For the homogeneous case, 
the ionization and recombination times corresponding to 
each cycle are identical and the largest cutoff (harmonic 
48a;) comes from the shortest pairs. 
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FIG. 3: (Color online) Dependence of harmonic order on the 
release time t' and the recombination time t of the electron for 
a nonhomogeneous field with e = 0.003 using the same param- 
eters as in Fig. (2] Panels (a) and (b) give the ionization and 
recombination times of SFA model, respectively, while panels 
(c) and (d) depict the ionization and recombination times of 
the classical calculations, respectively. The red dashed lines 
correspond to the harmonic cutoff. 

For the nonhomogeneous case, the general harmonic 
cutoff is extended in comparison to the homogeneous 
case, but the trajectories do not follow the same sym- 
metry as shown in Fig. [2l 

In here, the electron tunnels with two different canon- 
ical momenta given by Eq. (|3T|) : one leads to a higher 
cutoff, and the other to a lower one. For example, for 
the shortest pairs corresponding to the electron leaving 
at the field maxima the cutoff is at around harmonic 38w, 
while for the pairs corresponding to the electron leaving 
at the field minima it is at around harmonic 60a;. 

Furthermore, the cutoff become larger as we move from 
shorter pairs to the longer pairs. As it demonstrated in 
Fig. |4] the trajectories from the minimum and maxi- 
mum of the cycle moves towards each other until they 
collocate on each other at the field crossing. As e in- 
creases the solutions from positive and negative side of 
the cycle collapse faster to each other. It means there 
will be less available pairs of orbit for electron to pro- 
duce HHG. In case of the field with very large e even the 
long trajectory of the shortest orbits are not allowed as 
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FIG. 4: (Color online) Harmonic order as a function of the 
release times t' of the electron for the same parameters given 
in Fig. [S] but for pairs corresponding to the solutions from 
positive and negative sides of the first cycle. Panel (a) gives 
the SPA real part of the release times, while panel (b) depict 
the classical ionization times of these pairs. The dashed and 
solid lines correspond to the long and the short orbits, while 
the pairs with dot dashed lines do not have the well-known 
shorts and long pairs. 



we demonstrated it using the actual enhanced field of the 
nanostructure [32j . 



E. Spectra 

In this section, we compute HHG spectra using Eq. 

These results are depicted in Fig. [5] with red colored 
spectra presenting the homogeneous field and blue col- 
ored and black colored spectra showing the nonhomoge- 
neous case with e = 0.003 and e = 0.005, respectively. In 
all cases, the HHG cutoff is in good agreement with the 
trajectories analysis represented in the previous section 
and with the full 3D numerical calculations of Ref [l3 . 
These cutoffs are at around harmonics 45a;, 73a; and 93a; 
for e = 0, e = 0.003 and e = 0.005, respectively. More- 
over, it seems that both odd and even harmonics are 
present in the total spectra. At higher regime of the spec- 
tra, both odd and even harmonic have the same weight, 
while at the lower regime one sets of harmonics are a bit 
more dominant than the others. 



IV. CONCLUSIONS AND OUTLOOK 

We present theoretical studies of high-order harmonic 
generation in noble gases produced by spatially nonho- 




20 40 60 80 100 120 

Harmonic Order 



FIG. 5: (Color online) High-order harmonic spectra for hy- 
drogen atoms (Ip = 0.5 a.u.) interacting with a monochro- 
matic field of frequency uj = 0.057 a.u. and intensity / — 
3 X 10^*W/cm^. Red, blue and black colors depict the cases 
with £ = 0, e = 0.003 and e = 0.005, respectively. The HHG 
spectra obtained for e — and e = 0.005 are scaled for clarity. 



mogeneous fields. These fields are generated when differ- 
ent nanosystems, such as metal bow-tie nanostructures, 
metal and dielectric nanoparticles and nanotips, are ex- 
posed to a short and intense laser pulse. In this study, we 
have employed a numerical approach as well as performed 
an exhaustive theoretical analysis. 

For the numerical method, we used the numerical 
solution of the time dependent Schroedinger equation 
(TDSE) in reduced dimensions combined with the ac- 
tual functional form of the electric field extracted from 
finite element simulations to predict the harmonic spec- 
tra originated when Ar atoms interact with the plas- 
monic enhanced field in a metal bow-tie nanostructure. 
We observe an extension in the harmonic cutoff position 
that could lead to the production of XUV coherent laser 
sources and opening the avenue to the generation of at- 
tosecond pulses. For this particular nanostructure, this 
new feature is a consequence of the combination of a non- 
homogeneous electric field, which modifies substantially 
the electron trajectories, and the confinement of the elec- 
tron dynamics. 

In our theoretical analysis, we employed both classi- 
cal and semiclassical approaches to investigate the high- 
order harmonic spectra produced by spatially nonhomo- 
geneous fields of linear form. We show how the quan- 
tum orbits manifest themselves in these particular spa- 
tially inhomogeneous fields. In addition we predict that 
in nonhomogcneous fields, the electron tunnels with two 
different canonical momenta leading to a higher and lower 
cutoff. We also conclude that in the case of linear non- 
homogeneous fields, both odd and even harmonics are 
present in the HHG spectra. Within our model, we show 
that the HHG cutoff extends to the larger harmonics as a 
function of the inhomogcneity strength and this feature 
is consistent with the findings recognized using actual 
fields. 



8 



Acknowledgments 

We acknowledge the financial support of the MICINN 
project FIS2008-00784 TOQATA; ERG Advanced Grant 
QUAGATUA, Alexander von Humboldt Foundation and 



Hamburg Theory Prize (M. L.). This research has been 
partially supported by Fundacio Privada Gellex. We 
thanks Peter Hommelhoff and Samuel Markson for useful 
comments and suggestions. 



545 



163 



M. Protopapas, C. H. Keitel, and P. L. Knight, Rep. 
Prog. Phys. 60, 389 (1997). 

T. Brabec and F. Krausz, Rev. Mod. Phys. 72, 
(2000). 

M. Lein, J. Phys. B 43, R135 (2007). 
F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 
(2009). 

P. B. Corkum, Phys. Rev. Lett 71, 1994 (1993). 

M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L'Huillier, 

and P. B. Corkum, Phys. Rev. A 49, 2117 (1994). 

S. Kim, J. Jin, Y.-J. Kim, I.-Y. Park, Y. Kim, and S.-W. 

Kim, Nature 453, 757 (2008). 

I.-Y. Park, S. Kim, J. Choi, D.-H. L. Y.-J. Kim, M. F. 
Khng, M. I. Stockman, and S.-W. Kim, Nat. Phot. 5, 
677 (2011). 

P. Miihlschlegel, H.-J. Eisler, O. J. F. Martin, B. Hecht, 

and D. W. Pohl, Science 308, 1607 (2005). 

P. J. Schuck, D. P. Fromm, A. Sundaramurthy, G. S. 

Kino, and W. E. Moerner, Phys. Rev. Lett. 94, 017402 

(2005). 

A. Husakou, S.-J. Im, and J. Herrmann, Phys. Rev. A 
83, 043839 (2011). 

I. Yavuz, E. A. Bleda, Z. Ahun, and T. Topcu, Phys. 
Rev. A 85, 013416 (2012). 

M. F. Ciappina, J. Biegert, R. Quidant, and M. Lewen- 
stein, Phys. Rev. A 85, 033828 (2012). 
T. Shaaran, M. F. Ciappina, and M. Lewenstein, J. Mod. 
Opt. p. (in press) (2012). 

M. F. Ciappina, J. A. Perez-Hernandez, T. Shaaran, 
J. Biegert, R. Quidant, and M. Lewenstein, Phys. Rev. 
A 86, 023413 (2012). 

M. Sivis, M. Duwe, B. Abel, and C. Ropers, Nature 485, 
El (2012). 

S. Kim, J. Jin, Y.-J. Kim, I.-Y. Park, Y. Kim, and S.-W. 

Kim, Nature 485, E2 (2012). 

P. B. Corkum, (private communication) (2012). 

P. Hommelhoff, Y. Sortais, A. Aghajani-Talesh, and 

M. A. Kasevich, Phys. Rev. Lett. 96, 077401 (2006). 

M. Schenk, M. Kriiger, and P. Hommelhoff, Phys. Rev. 

Lett. 105, 257601 (2010). 

M. Kriiger, M. Schenk, and P. Hommelhoff, Nature 475, 



[22; 
[23; 
[24; 

[25; 

[26' 
[27 

[28; 

[29 



78 (2011). 

M. Kriiger, M. Schenk, M. Forster, and P. Hommelhoff, 
J. Phys. B 45, 074006 (2012). 

G. Herink, D. R. Solh, M. Guide, and C. Ropers, Nature 
483, 190 (2012). 

P. Saileres, A. L'Huillier, P. Antoine, and M. Lewenstein, 
Advances in Atomic, Molecular and Optical Physics, eds. 
B. Bederson and H. Walther 41, 83 (1999). 
A. L'Huillier and M. Lewenstein, Strong Field Laser 
Physics, ed. T. Brabec Springer Series in Optical Sci- 
ences (Springer, Berlin, 2008). 

Q. Su and J. H. Eberly Phys. Rev. A 44, 5997 (1991). 
J. L. Krause, K. J. Schafer, and K. C. Kulander, Phys. 
Rev. A 45, 4998 (1992). 

K. J. Schafer and K. C. Kulander, Phys. Rev. Lett. 78, 
638 (1997). 

S. S. Acimovic, Introduction to nanoparticle 
characterization in COMSOL (available from 



http://srdjancomsol.weebly.com 2011) 



[30] P. B. Johnson and R. W. Christy Phys. Rev. B 6, 4370 
(1972). 

[31] A. Thai, M. Hemmer, P. Bates, O. Chains, and J. Biegert, 
Opt. Lett. 36, 3918 (2011). 

[32] M. F. Ciappina, S. S. Acimovic, T. Shaaran, J. Biegert, 
R. Quidant, and M. Lewenstein, Opt. Express, (submit- 
ted), arXiv: 1204.591 Ivl (2012). 

[33] F. Siifimann and M. F. Kling, Proc. of SPIE 8096, 
80961C (2011). 

[34] F. Siifimann and M. F. Khng, Phys. Rev. B 84, 
121406(R) (2011). 

[35] C. H. Edwards and D. E. Penney, Differential Equations 
and Boundary Value Problems - Computing and Model- 
ing (Prentice Hall, 4th edition, 1999). 

[36] W. Gordon, Z. Phys 40, 117 (1926). 

[37] D. M. Volkov, Z. Phys 94, 250 (1935). 

[38] P. Salieres, B. Carre, L. L. Deroff, F. Grasbon, G. G. 
Paulus, H. Wahher, R. Kopold, W. Becker, D. B. Milo- 
sevic, A. Sanpera, et al.. Science 292, 902 (2001). 

[39] T. Shaaran, M. F. Ciappina, and M. Lewenstein, Phys. 
Rev. A 86, 023408 (2012). 



